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Abstract 


The paradigm of phenotypic plasticity indicates reversible relations of different 
cancer cell phenotypes, which extends the cellular hierarchy proposed by the classi¬ 
cal cancer stem cell (CSC) theory. Since it is still questionable if the phenotypic plas¬ 
ticity is a crucial improvement to the hierarchical model or just a minor extension to 
it, it is worthwhile to explore the dynamic behavior characterizing the reversible phe¬ 
notypic plasticity. In this study we compare the hierarchical model and the reversible 
model in predicting the cell-state dynamics observed in biological experiments. Our 
results show that the hierarchical model shows significant disadvantages over the re¬ 
versible model in describing both long-term stability (phenotypic equilibrium) and 
short-term transient dynamics (overshoot) of cancer cells. In a very specific case 
in which fhe fofal growfh of populafion due fo each cell fype is identical, fhe hier¬ 
archical model predicfs neifher phenotypic equilibrium nor overshoof, whereas fhe 
reversible model succeeds in predicfing both of them. Even though the performance 
of the hierarchical model can be improved by relaxing the specific assumpfion, ifs 
predicfion fo fhe phenotypic equilibrium sfrongly depends on a precondition fhaf may 
be unrealisfic in biological experimenfs, and if also fails fo capture the overshoot of 
CSCs. By comparison, it is more likely for the reversible model to correctly describe 
the stability of the phenotypic mixture and various types of overshoot behavior. 


1 Introduction 

The cancer stem cell theory has provided a hierarchical model of how diverse cancer cells 
being organized [[BEIEI . Similar to the stem cell theory in normal tissues, this hierarchical 
model assumes that a small number of stem-like cancer cells (termed cancer stem cells, 
CSCs) are capable of self-renewal and differentiation into other more committed cancer 
cells (termed non-stem cancer cells, NSCCs) but not vice versa. That is, CSCs are thought 
to be at the apex of this cellular hierarchy. However, some recent researches may extend 
this unidirectional relation of cancer cells. It has been reported that cancer cells can 
convert from NSCC phenotype to CSC phenotype (e.g. breast cancer [|4l[5l, melanoma 
DU, colon cancer DU and glioblastoma multiforme DU)- Furthermore, the interconversions 
among cancer cell phenotypes have also been found (breast cancer DU). All these works 
indicate reversible relations of different cancer cells (termed phenotypic plasticity, which 
has long been an issue concerned in bacterial populations ifTOl ). 

Very recently special attention has already been paid to the reversible cancer models 
with the phenotypic plasticity by theoreticians. In particular, dos Santos and da Silva ex¬ 
plained the variable frequencies of CSCs in tumors by establishing a model with stochastic 
cell plasticity DTTl |T2l . Leder et al investigated the model of reversible conversions be¬ 
tween the stem-like resistant cells (SLRCs) and the differentiated sensitive cells (DSCs) 
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Figure 1: Overshoot is a type of non-monotonie phenomenon [f20l that, starting from a 
state that is lower than the equilibrium level, the proeess first inereases above the final 
equilibrium level and then gradually deereases to it. 


in glioblastomas lfT3l . Wang et al showed how tumor heterogeneity arises in the model 
of eooperating CSC hierarehy with eell plastieity lfT4l . Zhou et al showed that the de¬ 
differentiation from NSCC phenotype to CSC phenotype was essential for explaining 
the transient inerease of the minority populations of CSCs observed in eaneer eell lines 
mm- Chen et al studied stoehastie models that eapture the transitions between en- 
doerine therapy responsive and resistant states of breast eaneer eells IfTTlI . Zhou et al 
investigated nonequilibrium dynamies with phenotype transitions of eaneer eells IfTSl . 
Jilkine and Gutenkunst studied the effeet of de-differentiation on time to mutation aequi- 
sition in eaneers m. 

However, it is still questionable if the phenotypie plastieity is a erueial improvement 
to the hierarehieal model or just a minor extension to it [121112^ . Thus a rigorous analysis 
on the eharaeteristies owned by the model with the phenotypie plastieity is neeessary for 
model validation. In this study, we try to investigate this issue by giving a eomparative 
study of the reversible model and the hierarehieal model. Note that Gupta et al studied 
SUM149 and SUM159 breast eaneer eell lines [|9l|, in whieh two interesting phenomena 
were observed: phenotypic equilibrium and overshoot. For the phenotypie equilibrium, 
they found that the breast eaneer eell lines will tend to a stable phenotypie mixture over 
time regardless of initial states. Similar results on the phenotypie equilibrium were also 
reported in ll23l [241 171. For the overshoot, they found that shortly after the eell sorting, the 
proportion of the minority subpopulation will inerease transiently above the final equi- 
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librium level, and then deerease to it (Fig. 1). Enlightened by these observations, we 
show that the reversible model is more eapable of eapturing the phenotypie equilibrium 
and overshoot than the hierarchieal model. Under a very speeific assumption that the 
total population growth due to each cell phenotype is identical, the hierarchical model 
can perform neither phenotypic equilibrium nor overshoot, whereas the reversible model 
succeeds in predicting both the phenotypic equilibrium and three types of overshoots 
(asynchronous, synchronous and oscillating overshoots). Even though the performance 
of the hierarchical model can be improved by relaxing the specific assumption, it is still 
not good enough to correctly capture the cell-state dynamics. On one hand, the pheno¬ 
typic equilibrium predicted by the hierarchical model strongly depends on the condition 
that the self-contributed growth rate by CSCs is faster than that of more committed cancer 
cells. However, it has been reported that this condition cannot be satisfied in some cancers 
[l25l[2^ . On the other hand, by fitting the two models to experimental data in [iHl, the hi¬ 
erarchical model can only fit the overshoot of NSCCs, but it cannot capture the overshoot 
of CSCs. By contrast, the reversible model can fit both of them. Therefore, our results 
generally imply that the reversible model shows distinct advantages over the hierarchical 
model in predicting both long-term and transient dynamics of cancer cells. 

The paper is organized as follows. The mathematical model is presented in Section 2. 
Main results are shown in Section 3, where we give a comparative study of the reversible 
and hierarchical models. Conclusions are presented in Section 4. 


2 Model 

In this section we describe the assumptions of our model. We model cancer as popula¬ 
tion dynamics of cancer cells, each cancer cell can be assigned to one of the following 
cell phenotypes: CSCq, NSCCi, NSCC 2 ,..., and NSCC„. Eet Xo{t), X„(f) be 

the cell numbers of CSCq, NSCCi, NSCC^ at time t respectively, then the model 
can be described as a cellular dynamical system of ^ = (Xo(f), Xi(f),..., X„(f))^. 
In particular, the model will reduce to a two-phenotypic model when n = 1. That is, 
besides CSC-phenotype, all the other cancer cells are grouped into one whole NSCC- 
phenotype. This CSC-NSCC model has extensively been investigated in previous liter¬ 
ature [12I1[II1[I21[I51[I1[I1. Here we pay more attention to the multi-phenotypic case 
(n > 2). As a starting point, the case with n = 2 should be a natural and favorable 
candidate for theoretical analysis. Meanwhile, motivated by Gupta et al @ where three 
phenotypes (stem-like, basal and luminal) were identified, we are concerned about CSCq- 
NSCC 1 -NSCC 2 model in light of both theoretical and experimental reasons. Unlike the 
three-phenotypic cell lineages investigated in |l27l [28l [23 where cancer stem cells, pro¬ 
genitor cells, and terminally differentiated cells are hierarchically cascaded, here we are 
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more interested in the eell population strueture investigated in [|9l|. In their work, CSCs 
ean differentiate into NSCCi and NSCC 2 , respeetively. They were not assumed to be eas- 
eaded. Moreover, by the phenotypie plastieity, they ean interconvert into each other and 
de-differentiate into CSCs. 

We now present the cellular processes included in our model. For CSCq, it can not 
only divide symmetrically into two identical CSCq daughters, but also divide asymmetri¬ 
cally into CSCo and NSCC^ (i = 1 or 2) That is, 

• CSCo ^ CSCo + CSCo. 

• CSCo ^ CSCo + NSC Cl. 

• CSCo ^ CSCo + NSCC2^. 


For NSCCi and NSCC 2 , besides the phenotypic plasticity, they can also perform symmet¬ 
ric divisions and cell death. That is. 


NSCCi 

01 

NSCCi + NSCCi 

NSCCi 

-02 

CSCo; 

NSCCi 

03 

NSCC2; 

NSCCi 

04 

0 ; 

NSCC2 

71 

■ NSCC2 + NSCC2 

NSCC2 

72 , 

CSCo; 

NSCC2 

73 , 

NSCCi; 

NSCC2 

74 

0 . 


Based on the above cellular processes, the dynamics of ^ = (Xo,Xi,X 2 )^ can be 
captured by the following ordinary differential equations (ODEs) 

dXo 


dt 

dXi 

dt 

dX2 

dt 


= ciiXo -|- 1^2X1 -|- ^2X2 
= 02-^0 + (/^l — 1^2 — — (d4)Ni + 73^2 

= 03-^0 + /dsNi -|- (71 — 72 — 73 — 74)-^2 


( 1 ) 


'it should be pointed out that, the symmetric differentiation “CSCo — >NSCCi + NSCCi” is also an 
important type of transition lf3^ that should be considered in the model. However, it can be shown that 
the main results we are interested in this work can still hold (with only minor changes in proofs) even if 
adding this transition into Eq. O- Therefore, to keep the model as minimal as possible, the symmetric 
differentiation is not included in presented model. 
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By letting 


72 

73 

7i - 72 - 73 - 74 


( 2 ) 


( ai /?2 

02 A - /^2 - /^3 - A 

«3 A 

Eq. Q can be expressed as 



(3) 


For this model, we have the following remarks: 

(1) Each element qij of Q can be seen as the per capita growth rate of phenotype 
contributed by phenotype. Let us take the first row of Q as an example. ai is the 
symmetric division rate of CSCq, thus the growth rate of Xq due to CSCq themselves can 
be expressed as aiXo; A and 72 are the de-differentiation rates from NSCCi and NSCC 2 
respectively, so the growth rate of Xq due to NSCCi and NSCC 2 is A-^i + 722 ^ 2 - In this 
way, it is easy to explain the biological meaning of “dXo/dt = aiX^ + A-^i + 72 -^ 2 ”, 
i.e. the rate of change of Xq is the sum of the growth rates contributed by all the three cell 
phenotypes in the population. Similarly, we can explain the other two equations in Eq. 

O. 

(2) Note that (^ 2 , 12 and 73 are the parameters associated with the phenotypic plasticity, 

the model will reduce to the hierarchical model by letting them be zero. Accordingly, Q 
will become to a lower-triangular matrix 


/ tti 0 ^ \ 

Q* = W= «2 A-A 0 . (4) 

\ 03 0 7i - 74 / 

{3) Let W = Xq -f Xi -f X 2 be the total number of the population, then 
dW 

—— = (ai -f 0:2 + OiQ)XQ + (A — A) 2 fi + (71 ~ 74)2^2- ( 5 ) 

dt 

Note that the coefficient of each Xi on right-hand side of Eq. Q is just the the sum of the 
corresponding column in Q, i.e. the per capita growth rate of the whole population due to 
Xi. In particular, when 


tti -f 02 + as = A - A = 7i - 74 = (6) 

the growth rate of the whole population due to each phenotype will be identical. This 
special case implies meaningful biological significance and has been reported in some 
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cancer cell lines (breast eancer Q and eolon eaneer IfTSl l. In this ease Eq. ([5]) will reduee 
to a very simple form 


dW 

dt 




that is, W will grow exponentially with constant rate. We will show that, even in this spe- 
eial case the reversible model and the hierarchieal models are quite different in predieting 
the phenotypic equilibrium and overshoot. For eonvenienee, we denote 


Ai — Oil -\- 0.2 0 ^ 3 ) 


(7) 


A2 — fdi — jdi, 


( 8 ) 


^3 = 7i - 74- (9) 

(4) Note that all the parameters are assumed to be eonstant. In other words, there is no 
feedback control included in presented model. This is of eourse a simplifieation of the 
biologieal faets. To model real biologieal systems, feedbaek meehanisms are indispens¬ 
able lf34l [35l [36l |28l . However, we will show that our model is not overly simplistie to 
eharacterize the phenotypie plastieity in eomparison to the eellular hierarchy. 

In the remaining of this section, we will convert the cell number equations Eq. ([T]) 
into the following cell proportion equations Eq. (fTOl) . This is because in reality, relative 
eell numbers, rather than absolute eell numbers, are usually measured by fluoreseenee- 
aetivated eell sorting (FACS) experiments. Eet xq = Xq/W, Xi = Xi/W and X 2 = 
X 2 /W, by variable substitutions in Eq. ([T]) we have 


" dxo 
dt 

dxi 

—- 

dt 
dx2 
' dt 

Note that Xq 
1 - Xo - Xi, 

{ dxQ 

dt 

dxi 

IF 


— aiXf) + f 32 Xi + 722^2 ~ Xo{AiXq + A2X1 + ^43X2) 

= a2Xo + {/3i - 132- ^3- Mxi + 733^2 - a:i(AiXo + A2X1 -f A3X2) (10) 

= 0:3X0 (dsXi + (71 — 72 ~ 73 ~ 74)2^2 ~ X2{AiXo + A2X1 + A3X2) 

+ Xi -1- X 2 = 1, the freedom of Eq. (fTOl) is actually two. By letting X 2 = 
Eq. (fT^ ean be transformed into 

= -(Ai - A3 )xo - (A2 - A 3 )xoXi + A4X0 + A5X1 + 72 

(11) 

= -(A2 - A3)xi - (Ai - A3 )xoXi + AqXo + A7X1 + 73 
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where 


^4 = ai - (71 - 74 ) - 72 , 


( 12 ) 


A^ = (^2- 72, 


(13) 


Aq — a 2 — 73 , 


(14) 


Ai = (/3i - /32 - /^s - /34) - (71 - 74 ) - 73- 


(15) 


When /92 = /^s = 72 = 73 = 0, /.e. the phenotypie eonversion rates are set to be zero, we 
have 


— —{Ai — A^)xl — {A2 — A^)xqXi + ^4X0 + ^53^1 


dxo 
dt 

= —{A2 — ^3)3^1 ~ (^1 ~ ^ 3 )^ 02 ^! + A ' qXq + A'rjXl 


(16) 


where 


A '^ = ai - (71-74), 


(17) 


^5 = 0, (18) 

4 = a2, (19) 

^7 = (/5i-/34)-(7i-74). ( 20 ) 

We term Eq. (fTTl) the reversible model, and term Eq. (fT^ the hierarchical model. The 

eomparison of these two models is the main task in the following section. 

In this section we try to give a comparative study of Eqs. (fTTl) and (fUf). 

One trivial difference between the two models happens when there is no CSC in the 
initial population. In this case, since in the hierarchical model CSCs cannot be sponta¬ 
neously created by other cell phenotypes, there is no CSC in the population all the time, 
whereas new-born CSCs can be generated from NSCCs in the reversible model. However, 
in reality it is quite impossible to completely purify cancer cell lines without any CSCs 
by cell sorting. Based on this fact, we only discuss the cases with positive initial states 
of CSCs afterwards. Both long-term stable behavior and short-term transient dynamics 
are investigated when comparing the two models. Enlightened by S, we are particularly 
concerned about their predictions to the phenotypic equilibrium and overshoot. 
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2.1 Special case 


Before exploring the general models in Eqs. (fTTI) and (fT^ . we first diseuss a very speeifie 
case with ^41 = ^42 = A^. Investigating this special case can provide some valuable 
insights to the general analysis. 

For the hierarchical model in this case, 


{ dxp 

dt 

dxi 

Ht 


A 4 X 0 


-46X0 


( 21 ) 


It is easy to give their explicit solutions 

f Xo(t) = Xo{0)e~A^~^°‘3)t 

[ Mt) = ^-o(0)(l - +ii(0) 


where xo(0) and xi(0) are the initial states. So we have 


1. Since xo{t) and xi{t) are both monotonic functions of t, they cannot perform over¬ 
shoot. 


2. Note that 


lim Xoit) = 0, 

i —^~|“00 


0(9 

lim xi{t) = -^-Xo(0) +xi(0), 

t ^+00 0!2 C^S 

CSCs proportion will eventually be zero, and the equilibrium proportions of xi 
depends on the initial states. Note that the phenotypic equilibrium corresponds to 
the stabilization of the phenotypic mixture that is independent of the initial states 
[|9l , the hierarchical model fails to predict the phenotypic equilibrium. 


We turn our attention to the reversible model 


dxo . . 

—7— — A4X0 + ^5X1 + 72 
dt 


dxi 

dt 


= AqXo + A-iXi -f 73 


(23) 


For the phenotypic equilibrium, we have the following theorem: 
Theorem 1. There exists unique stable fixed point E = (xg, x^) in Eq. 

* ^ -4573 - -4772 ^ ^ -4672 - -4473 

° AA-- 45-46 ’ ' 7I42I7 - 2I52I6 


m, where 

0 . 
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Table 1: The overshoot predieted by the speeial reversible model 



Parameters 

Behavior 

1 

A > 0 

no overshoot, asynchronous and synchronous overshoots 

2 

A = 0 and A is diagonalizable 

no overshoot 

3 

A = 0 and A is not diagonalizable 

asynchronous and synchronous overshoots 

4 

A < 0 

oscillating overshoot 


The proof is put in El Theorem [Hindioates that the phenotypie equilibrium naturally 
holds in the speeial reversible model. We now eonsider whether this model ean perform 
overshoot. The main results are listed in Table [T| (see |B] for more details). In partieular, 
Fig. 2 illustrates the predietions by ease 1 (A > 0): 

• Non-overshoot (leftpanel). Both xq and xi are monotonie funetions. 

• Asynchronous overshoot (middle panel). Either xq or xi ean perform overshoot 
(The figure shows an example of a;o-overshoot), but they eannot perform overshoot 
simultaneously. 

• Synchronous overshoot (right panel), xq and xi ean perform overshoot simultane¬ 
ously. 

Case 2 (A = 0 and A is diagonalizable) is the star ease (see E]), overshoot ean never 
happen. For ease 3 (A = 0 and A is not diagonalizable), it belongs to the nodal ease (see 
|B]), both asynehronous and synehronous overshoots ean happen. Besides, in ease 4 (A < 
0 ), the model performs damped oseillatory dynamies, eorresponding to an interesting 
oscillating overshoot (Fig. 3). 

Aeeording to the above eomparisons between Eqs. (I2TI) and (1231) . the reversible model 
shows rieher dynamics than the hierarchical model in performing both the phenotypic 
equilibrium and overshoot. 

2.2 General case 

We now remove the constraint of Ai = A 2 = A 3 and further consider the general case in 
this section. For the general hierarchical model Eq. (fT^ . the following theorem states the 
condition under which the phenotypic mixture will finally stabilize: 

Theorem 2. There exist three singular points in Eq. Eq = (0, 0), E[ = (0, 1) and 
E 2 = (xgja;*), where 

* ^_ ^4^6 _ * ^ (A4 — Aj)xl 

^ (ai-a',)(a;-a') + a'A'’ ° A' • 
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Non-overshoot 


Asynchronous overshoot 


Synchronous overshoot 




Figure 2: Three types of predictions by the special reversible model in case 1 (A > 0). 



Figure 3: Oscillating overshoot: Compared to the simple uphill-downhill overshoot in Fig 
1, oscillating overshoot shows multiple uphill-downhill movements with damped ampli¬ 
tude. 
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Assume that ai > /3i — [5^ and ai > 71 — 74, Ei^ = is the unique stable fixed 

point. 

The proof is put inO Note that ai, fii — and 71 — 74 are the diagonal elements of 
the matrix Q* in Eq. dH), representing the per eapita growth rates due to the eorresponding 
cell phenotypes themselves. Theorem [21 thus indicates that the general hierarchical model 
is capable of predicting the phenotypic equilibrium only provided that the self-contributed 
growth rate by CSCs is larger than that of more committed cancer cells. However, this 
condition may not be satisfied in reality. It has been reported that more committed cancer 
cells can have faster cycling time than stem-like cancer cells in some cancers [[25l [2^ 
IT4]| . In contrast, the following theorem shows the phenotypic equilibrium of the general 
reversible model Eq. (ITOl) : 

Theorem 3. The reversible model Eq. (|70|) has unique positive stable fixed point^ 

The proof is put injDl Theorem |3] shows that the result in the special case (Theorem [T]) 
can be extended to the general case. By comparison of Theorems |2] and IH we can see that 
it is more likely for the reversible model to correctly capture the phenotypic equilibrium. 

In the remaining of this section, we discuss the overshoot performances of the two 
models by fitting them to the cell-state dynamics in [|9l. We fitted Eqs. (fT^ and (fTOl) to 
the SUM159 data shown in the Eig. 3 of [|9|. Erom Eig. 4, we can see that the reversible 
model fits the data better than the hierarchical model. In particular, note that both the 
overshoots of luminal cells and stem-like cells were presented in the SUM 159 data (see 
the left panel and right panel in the Eig. 3 of ||9l respectively). Our result shows that, the 
reversible model can capture both of them, whereas the hierarchical model can only fit the 
overshoot of luminal cells but it fails to capture the overshoot of stem-like cells. In the left 
panel of Eig. 4, both the blue solid line (hierarchical model) and blue discrete triangle- 
marker line (reversible model) fit the overshoot of luminal cells. Meanwhile, in the right 
panel of Eig. 4, the red solid line (hierarchical model) shows a monotonic way up to 
the final equilibrium (non-overshoot), but the red discrete triangle-marker line (reversible 
model) predict the overshoot of stem-like cells. 

Actually, it has already been reported in previous literature that the overshoot of 
NSCCs proportion can happen in classical CSC model [[371l . but the overshoot of CSCs 
proportion we think has added significance. Note that in the hierarchical model the in¬ 
crease of CSCs relies only on their self-renewals. When the initial fraction of CSCs is 
very limited, the transient increase of CSCs proportion should be very slow (constraint by 
the limitation of cell division cycle), so intuitively it is quite unlikely for the hierarchical 
model to capture the overshoot of CSCs proportion. In contrast, the de-differentiation 

^We should point out that, to complete the hnal proof, technically it is necessary to add a small pertur¬ 
bation to the initial state in few cases, seelDl 
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Figure 4: Predictions of the general hierarchical model and reversible model. We used 
least square method to fit the two models to the SUM 159 data in [0. The left panel shows 
the purified stem-like cell case, where the initial states of stem-like cells, basal cells and 
luminal cells were assumed to be 99%, 0.5% and 0.5% respectively. The right panel shows 
the purified luminal cell case, where the initial states of stem-like cells, basal cells and 
luminal cells were assumed to be 0.5%, 0.5% and 99% respectively. In both panels, the 
solid lines are the predictions by the hierarchical model, and the discrete triangle-marker 
lines are the predictions by the reversible model. 


from NSCCs can effectively speed up the accumulation of CSCs, it is more likely for 
the reversible model to capture the overshoot of CSCs proportion. In other words, the 
overshoot of CSCs should be a better candidate than that of NSCCs to characterize the 
phenotypic plasticity. Our numerical result just supports this idea. 


3 Conclusions 

In this paper, we have tried to characterize the phenotypic plasticity by giving a compara¬ 
tive study between the reversible model and hierarchical model. According to our results, 
the reversible model shows richer dynamics and better predictions to experimental phe¬ 
nomena than the hierarchical model. 

It should be noted that, to focus our attention to the reversible phenotypic plasticity, 
the presented model does not include the biologically complex mechanisms such as feed¬ 
back controls and time delays. It is undeniable that, by including the feedback and delay 
mechanisms, the hierarchical model can also show complex dynamics that are beyond the 
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reach of the presented model [lMl[39l|40l[34l[^. For example, Bernard et al showed that 
oscillations can occur in their hierarchical model with feedbacks and delays (SSl- Stiehl 
et al showed that the hierarchical model of hematopoiesis with nonlinear feedback can 
reproduce the overshoot and damped oscillations observed in clinical data [[39l. We also 
think that the interplay between the feedback (or delay) mechanisms and the phenotypic 
plasticity should be a very interesting issue. As a staring point, it is suggested by our 
presented model that, the phenotypic plasticity facilitates the heterogeneity of cancer in 
two ways: it helps to sustain the coexistence of multiple phenotypes in cancer, and it 
also accelerates the recovery process of cancer stem cells in the population. These find¬ 
ings may shed some lights on the researches of the models incorporating the biological 
complexities with the phenotypic plasticity. 

Moreover, the presented model is deterministic, whereby stochastic effects are not 
accounted for. Therefore, the branching model iHTI concerning both the stochastic phe¬ 
notypic conversions and proliferations of cancer cells is another interesting research di¬ 
rection. Finally, it is worth noting that, in addition to the long-term stabilization, the 
short-term transient dynamics also deserve special attention in the study of biological 
population dynamics W2\ . In this work we have seen the significance of the overshoot in 
characterizing the phenotypic plasticity. It is interesting to explore more types of transient 
dynamics and their biological significance in future researches. 
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A Proof of Theorem U 


Theorem (pagelH TheoremU]). There exists unique stable fixed point E = (xg, xj;) in Eq. 
t l23l) . where 


^ _ A573 — A772 ^ _ ^672 — A473 

° ~ A4A7 - A5A6 ’ ~ A4A7 - A5A 


Proof. Let 


f A4X0 + A5X1 -f 72 — 0 

\ AgXo + A7X1 -f 73 = 0, 


(24) 
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where A 4 — — («2 + Q^ 3 + 72 ) < 0, ^7 — — (/32+/?3+73) < 0, ^5 — /32~72> ^6 — Q!2 — 73. 

Note that 


A4AY — A^Aq > A573 — A772 > 0, ^44^47 — A^Aq > Aq'^2 — ^473 > 0, 


it is easy to see that E = is the only fixed point. To proof the stability of E, let 

xq = Wo + Xq, xi = wi + xl, then Eq. (|2^ is transformed to 


dwo 

dt 

dwi 

dt 


= A4W0 + A^wi 
= AqWq + AjWi 


Let w = {wq, Wi)'^ and 

A={ 

\ Aq Aj 

Eq. (I 25 ]) ean be expressed as 


dw 

dt 


= Aw. 


E' = (0, 0) beeomes the only fixed point and its eharaeteristie equation is 
— {A 4 + Aj^X + A 4 A 7 — A^Aq = 0 . 


(25) 


(26) 


(27) 


Sinee 

Ai + A 2 = A 4 + ^47 < 0 , A 1 A 2 = A 4 AY — A^Aq > 0 , 

the real parts of Ai and A 2 are both negative. Henee E' is stable in Eq. (l25l) . and then E is 
the unique stable fixed point of Eq. (|2^ . The proof is eompleted. □ 

B The overshoot predicted by the special reversible model 


Let 


A — {A 4 — A-jY + dAsAg (28) 

be the diseriminant of (1271) . Sinee it has been shown that the real parts of Ai and A 2 are 
negative, the speeial reversible model Eq. (|2^ ean be elassified into the following four 
eases (Eor the terminologies used below, please refer to Chapter 5 in [|43l): 
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Figure 5: The phase diagrammatie sketeh of the ease that A = 0 and A is diagonalizable. 




Figure 6: The phase diagrammatie sketeh of the ease that A = 0 and A is not diagonal¬ 
izable. (a) Asynehronous overshoot: When starting from Pi, xq ean first exeeed Xq and 
then deerease to it (overshoot), whereas xi eannot perform overshoot. However, when 
starting from P 2 , Xi ean perform overshoot, but Xq eannot perform overshoot. That is, 
either xq or xi ean perform overshoot, but they eannot perform overshoot simultaneously, 
(b) Synehronous overshoot: When starting from Pi, both xq and xi ean exeeed E and 
then deerease to it. That is, they ean perform overshoot simultaneously. 


1. A > 0. In this ease the matrix A in Eq. (|2^ has two unequal negative eigenvalues. 
Thus E' is the nodal point of Eq. (|25]) . eorrespondingly E is the nodal point of Eq. 
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(|2^ . The solution can be expressed as 


(29) 


where 


P 


Pll Pl2 
P2I P22 


(30) 


is a matrix determined by the coefficients of Eq. (l25l) . both ki and ^2 are determined 
by the initial states. Different types of transient dynamics can happen for different 
parameter ranges. By numerical simulations, it is shown in Fig. 2 that the model 
can perform non-overshoot, asynchronous overshoot and synchronous overshoot 
respectively. 

2. A = 0 and A is diagonalizable. Since E is the star point in this case (Fig. B.5). 
There is no overshoot of xq and xi. 

3. A = 0 and A is not diagonalizable. In this case E is also called nodal point. 
From Fig. B. 6 , we can see that the model can perform both asynchronous and syn¬ 
chronous overshoots. If (A 4 — A-j)jA^ > 0, the model can perform asynchronous 
overshoot (Fig. B. 6 (a)). If {A 4 — A'j^jA^ < 0, the model can perform synchronous 
overshoot (Fig. B. 6 (b)). 

4. A < 0. A has two conjugate complex eigenvalues with negative real part. E is the 
focal point, resulting in oscillating overshoot (Fig. 3). 


C Proof of Theorem |2| 

Theorem (page [l2l Theorem [2]). There exist three singular points in Eq. (E^).- Eq = 
(0, 0), E[ = (0,1) and E 2 = {xq, Xi), where 



A' A' 

ZI4/I6 



Assume that ai > (3i — and ai > 71 — 74 , E '2 = {xq,xI) is the unique stable fixed 
point. 


Proof. To obtain the fixed points of Eq. (fT^ . we consider 

f -(A( - A'^)xl - {A '2 - A'fjXoXi -f A'^xo = 0 
\ -(A 2 - A'f)x\ - (A( - A'fixoXi -f A'qXq + A'rjXi = 0 


,2 

0 

,2 


(31) 
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For the first equation — {A'^ — ^3)0:0 — [A'^ — A'^)xqXi + A'^Xq = 0 , we have either Xq = 0 
or — {A[ - A'^)xq — {A2 - A'^)xi + ^4 = 0. 

(i) When xq = 0 , based on the seeond equation we have xi = 0 or xi = A'-j/{A'2 — A'^) = 
1, that is, we have two fixed points 

E' =(0,0), = (0,1). 

(ii) When —{A'^ — A'^)xq — {A'2 — A'.^)xi + A4 = 0 , the fixed point is 

-^2 = (%) ^ 1 )) 


where 


Xi = 


A' A' 


Xcx = 


(^4 ^ 7 ) 3^1 


Ai 


[A'^ A'-j)[A'^ A'^) + A'qA'^ 

Sinee ai> (3i — 134 and ai > 71 — 74 , we have 0 < Xq < 1 and 0 < X 4 < 1. Let 

/(xo,Xi) = -(A'l - A'^)xl - {A2 - A3)xoXi + A4X0, 

g{xo, xi) = -{A 2 - A'^)xl - {A[ - A3)xoXi + AgXo + A'jXi, 

then 

/xo = “^(A^ - A'^)Xo - {A2 - A'^)x 4 + ^4, = -(A2 - A'^)Xo, 

9 x 0 ~ ~(^1 ~ ^3)^1 + 9 x 1 ~ ~ ^3)^1 ~ (^1 ~ ^3)^0 + Ay 

The Jacobian matrix of Eq. (fT^ is given by 


f' f 

J — I JXQ JXl 

\9xo 9x1 


Let us first consider the stability of Eq = (0,0), 


= a 4 ; 


its characteristic equation is 


(A-Al)(A-4) = 0. 


Note that 


ai> / 3 i- (34, tti > 7i - 74, 

then Ai = ^4 > 0, implying that Eq = (0, 0) is unstable. 
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Similarly, for E[ = (0,1) 


J{E[) 


( A',-A', 0 \ 

U-^'i + ^3 -A',) 


Xi=A'^- /1'7 > 0, E[ 
For 


where 


(A — (^4 — A7))(A + A'^) — 0, 
(0,1) is also unstable. 





fxi\ 

9 x 1 J ’ 


fxo = fx 




iJy 


/xi = fx 


XI ^ 


iJy 


9x0 = 9'xoi^o^f)^ 


9x^ = 9'xAXo,x] 


Then we have 

X ifxQ T 9xfX T fxo9xi fxi9xo 0, 

A = [—(Ai + A2)]^ — 4 A 1 A 2 = A'^ > 0, 

Ai + A 2 = A'y — 2^4 < 0, 

A 1 A 2 = ^ 4(^4 — A'rj) > 0. 

Then 

Ai < 0 , A 2 < 0 . 

It is shown that E 2 is a stable nodal point, and then it is the unique stable fixed point in 
Eq. (fT^ . The proof is eompleted. □ 


D Proof of Theorem |3| 

Theorem (paged^ Theorem [3]). The reversible model Eq. HUi has unique positive stable 
fixed point. 

Proof. Our proof is valid for general multi-phenotypie eases, so in this seetion the dimen¬ 
sion of the matrix Q in Eq dH) is not restrieted to three any more, but n in general. Without 
loss of generality, we only prove the theorem when Q is positive (all the elements of Q are 
positive)!^. Aeeording to the well-known Perron-Frobenius theory (see ehapter 1 in EH), 

^Note that the off-diagonal elements of Q are positive, Q + tI (I is the identity matrix) can be positive 
provided t is large enough. It is easy to show that when A is an eigenvalue of Q -f tI, A — r is an eigenvalue 
of Q accordingly. That is, the eigenvalue structure of Q is the same as that of Q + rl. So we just need to 
consider Q + tI instead of Q. 
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Q has a real eigenvalue Ai satisfying ReX < Ai for any other eigenvalue A of Q (called 
the Perron-Frobenius eigenvalue) and Ai is simple (a simple root of the characteristic 
equation of Q). 

Since Ai is simple, the solution of Eq. ([T]) can be expressed as 


f, = 


j=2 1=1 


'hi 


i=l 




(32) 


where Ai, A 2 , • • • Am are the different eigenvalues of Q, mj is the algebraic multiplicity 

of Xj, u is the normalized {ui + U 2 + ... +«„ = !) right eigenvector of Ai, r/. is the 
corresponding eigenvector of Xj, Cjj is determined by initial states. Suppose Ci^i 7 ^ 0, 
since ReAj < Xi {i ^ 1), 


3^ 


7=2 1=1 i=i 


{X,+ 

Noticing that Perron-Frobenius theory ensures that all the components of u are positive, 
i.e. u > 0. 

Before we complete the proof, we need to discuss ci^i = 0. In this case, the above 
method does not work. However, since fluctuations are inevitable in real world, ci 1 = 0 
will hardly happen in reality. To show this, let f = 0 in Eq. (l3^ 

n n-j ^ 

Cl,lM + XI 

i =2 1=1 

This is a linear equation of Cjj. For ci^i we have 

_ det|5*| 

“ det|5| ’ 

where B = [u rf^ ^], B* is just B with its first column replaced by (^)^. It 

is easy to add a small perturbation ev to (Xq)'^, so that all the columns of B* are linear 
independent, hence ci,i 7 ^ 0 . 

□ 


_^_ u _ 

... + Xn) / Ui U2 -\- ... + Un 




^ 3^ 

"* = IE = 
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